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DETECTING  A  POINT  RADI/ TOR  IN  THE  PRESENCE  OF 
NONGAUSS  IAN  (INTERFERENCE 


[Reznik,  M,  Kh.,  Obnaruzheniye  T.-chechnogo  Izlychatelya  v  Prisutstvii 
Fonovykh  Pomekh  Negaussovogo  Tip.t,  Optiko-mekhanicheskaya  promyshlen- 
nost',  No.  3,  1972,  pp.  3-6,  Rus  ian] 


An  optimal  method  of  detecting  a  pulse  signal  in  the  [3 

presence  of  background  interference  of  the  kind  presented 
by  cloudy  sky  radiation  characterized  by  a  nongaussian 
distribution  law  is  described. 


Many  problems  in  optical  lo<  ation,  navigation  direction  finding, 
and  the  like,  are  connected  with  the  detection  of  .  signal  from  a  point 
radiator  against  a  background  of  natural  formation: . 

In  works  known  to  date  [1-3]  the  effect  of  ba«  kground  interference 
was  taken  into  account  within  the  framework  of  cor  elation  analysis, 
leading  to  linear  methods  of  analyzing  the  electri  al  signal  arising 
in  the  radiation  receiver.  Experience  in  the  use  f  linear  systems 
shows  that  in  the  case  where  signals  are  received  rom  point  objects, 
t.e.,  signals  having  the  form  of  Dirac  &  -functio  s,  linear  systems 
are  least  protected  from  the  effects  of  interferes  e  whose  values  change 
spasmodically. 

The  construction  of  a  statistically  optimal  s:  stem  consists  first 
of  all  in  the  selection  of  a  background  model  that  reflects  its  spas¬ 
modic  structure.  The  simplest  such  model  is  a  bad  ground  model  in  the 
form  of  a  population  of  spots  (clouds)  randomly  di»  tributed  against  a 
"clear  sky"  and  having  random  (but  constant  within  the  limits  of  each 
spot)  brightness  values.  ~We  shall  assume  that  in  any  cross  section  of 
such  a  background  the  distribution  of  boundaries  (transitions  from  dark 
to  light  and  back  again)  is  similar,  i.e.,  that  the  background  is  homo¬ 
geneous  and  isotropic,  and  that  the  boundary  distribution  is  Poissonian, 
i.e. , 

P(n.  r)  I- (>r)« exp (->«),  (1) 

where  p(n,r)  is  the  probability  that  in  the  intervj 1  r  w  11  be  equal 
to  n  Intersections;  v  is  the  average  frequency  of  intersections  in  a 
single  segment. 

It  is  useful  to  distinguish  this  model  from  the  anisotropic  two- 
level  background  model  in  the  form  of  a  "chess  board,"  examined  in  works 

(3-41. 
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We  shall  first  examine  the  interference  signals  that  form  during 
search  scanning  (in  the  case  of  uniform  displacement  of  the  radiation 
receiver  along  the  x  axis)  when  a  background  inhomogeneity  boundary  [4 
is  intersected  in  the  assumption  of  ideal  optics.  The  shortcomings 
of  optics,  or  aberration  distortions,  we  shall  consider  later. 

We  shall  assume  that  the  radiation  receiver  is  in  the  form  of  a 
Hxh  (H>h)  rectangle  whose  angular  dimensions  are  much  less  than  the 
mean  angular  dimensions  of  the  background  inhomogeneity,  i.e., 


v//  i.  (2* 

Granting  this,  the  inhomogeneity  boundary  can  be  prese  ted  in  the  form 
of  a  direct  boundary  with  inclination  0  to  the  scannin  direction.  The 
signal  received  from  the  interference  $  (x),  received  as  the  result  of 
boundary  intersections,  for  a  considerable  part  of  its  change  Ls  des¬ 
cribed  as  a  linear  function  (as  shown  in  figure  1,  wit  recei\  ir  dis¬ 
placed  from  position  1  to  position  2).  This  makes  it  ossiblc  to  ap¬ 
proximate  its  increase  or  decrease  with  a  linear  funct  on  thrc  lghout 
the  entire  sector  of  change.  The  signal  from  the  poin  object  has  the 
for-n  of  a  rectangular  pulse  having  width  h. 

We  shall  eliminate  the  constant  comporent  present  in  the  Interfer¬ 
ence  either  by  means  of  linear  transformation  or  by  meins  of  <■  second 
receiver  located  beside  the  first  and  connected  in  an  <>pposit<  way.  There¬ 
fore  the  resulting  signal  from  the  interference  «l»,  (jc)  -?(.v)  ?(.v-  -//l 

is  realized  in  a  form  consisting  of  nonintersecting  isosceles  trape¬ 
zoids  located  along  the  abscissa  (capable  of  degenerating  inti  tri¬ 
angles)  of  positive  and  negative  polarity,  random  extension  ai  i  ampli¬ 
tude  (but  limited  on  top  along  the  modulus,  since  maximal  bad  ground 
brightness  drops  are  limited)  and  randomly  separate  from  each  >ther. 

It  is  important  to  note  that  the  frontal  lengths  of  th  i  trape.  iidal 
formations  are  precisely  equal  to  the  width  of  the  rectangula  pulse. 

The  signal  from  the  point  object  is  obtained  in  t  le  form  of  rec¬ 
tangular  pulses  of  different  polarity  having  wic'th  h  a  jacent  in  the 
base. 


The  effect  of  aberration  distortions  manifests  itself  in  the  "smooth¬ 
ing"  of  the  processes  obtained.  As  a  result  of  this  "smoothii  the  rec¬ 
tangular  pulse  assumes  a  bell-shaped  form.  The  function  deso  Ibing  its 
change  as  a  function  of  the  time  argument  is  designates)  s  (t)  ind  it  is 
assumed  that  at  t,  >  A/2  in  (/)  0  Aberration  distortions  o?  opi  .cs  can  be 

treated  as  a  kind  of  linear  transformation,  consequently,  the  "smoothed" 
interference  fronts  will  as  before  have  a  length  equa?  to  that  of  signal 
S'  and  be  described  by  functions  that  are  integrals  of  so(t). 
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Figure  Is  Change  in  receiver  position  when  inter-  acting 
inhomogeneity  boundary:  a)  0<8*,  b)  0  >  0*. 

The  random  interference  time  function  $  (t)  can  be  presented  in 
the  *001  4  (t)  •  F(t)A(t) ,  where  function  F(t),  normed  to  absolute 

maxi-ium  value,  is  the  same  type  of  function  as  <|>(t),  while  A(t)  >  0 
is  a  spasmodic  process  with  a  uniform  distribution  of  values  :ln  segment 
[o,Ajiax]  that  changes  its  values  only  at  those  time  moments  when  process 

F(t)  deviates  from  the  zero  state.  For  ease  of  later  computations,  we 
shal.  consider  that  process  A(t)  assumes  only  a  finite  number  of  values 
0  <  fli  <  Qj  <  . . .  <  ar  *  .  />  '$>  1  •  is  a  whole  number  and 


P(A  (t)  -  Oj)  —  ~ ;  /=■•!,  2 - p.  (.'■!> 

One  of  the  possible  realizations  of  interference  ^(t)  :ogether 
with  functions  A(t),  F(t),  and  sQ(t)  is  depicted  in  figure  2. 

The  signal  from  the  point  object  will  become  bipolar  and  will  be 
determined  by  function  fcs(t),  where  fc  is  the  signal  amplitude  unknown 
to  the  observer. 

Radiation  receiver  noise  is  considered  white  and  normal. 
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Since  we  are  assuming  that  the  radiator,  subject  to  detection,  is 
a  point  radiator,  signals  from  the  target  and  background  will  be  addi¬ 
tive  and  the  total  signal  will  be  determined  as 

^  (/)/'(/)-!•«  (0  4  :*(/). 

if  there  is  a  signal,  and 
A  (t)!(t)+n(t), 
if  there  is  not  a  signal. 
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Figure  2:  Realization  of  random  functions 
$(t),  F (t) ,  and  A(t). 


Without  disturbing  the  generality  of  the  examinat  on,  we  solve  the 
problem  in  the  discrete  form,  assuming  that  the  observ  tion  of  process 
y(t)  occurs  at  discrete  time  moments 

-  f.v;  /,  =  /,  (/-  l)  Sf 

and 

y«Ky(W 

are  values  if  the  observed  process  at  reterence  point* .  We  shall  select 
the  separation  interval  so  that  there  is  a  sufficient  number  of  readings 
for  the  duration  of  the  signal,  i.e.,  so  that  b/  t  »  k>  1.  Receiver 
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noise  will  be  given  by  the  sequence  nj  of  normal,  mutually  independert 
values  with  parameters  (0, <T  ) .  In  examining  the  problem  in  a  time- 
discrete  sense,  and  assuming  that  k  %  1,  it  is  natural  to  assume  that 
-  both  the  signal  arrival  time  and  the  time  moments  of  transition  of  pro¬ 

cess  F(f)  from  one  set  of  states  to  another  correspond  with  selected  t^ 
readings.  Whence  it  follows  that  interference  values  F(t)  are  completely 
i  determined  by  the  set  of  4k-l  of  f.  numbers,  where  i  «  -(2k  -  1)  ... 

(2k  -  1),  which  are  computed  with  the  formulas 


/o  =  0, 

-  ip+ut  i-y 

*•  (y)dyi 

i=l,  2  ...  k, 

/=/»-/;  i~k+  1,  k  +  2  ...  2*-l;  (5) 

fi~  ~f-i,  i<0. 

The  v  ilues 

/,.  /,-••/*  U  !./->•••  /-  <*-»>  and 

/  .  t.  /tJ-j  .  .  .  fik-l  (/-I ,  /-  ?  •  •  •  /-  <»”•>) 


//=  J 


-C/2 


correspond  respo  rtively  to  increasing  (falling  in  the  negative  region) 
and  falling  (increasing  in  the  negative  region)  interference  pulse 
fronts;  f^  and  f  ^  are  extremal  internal  values. 

The  s  .gnal  -/ill  be  given  by  vector 

St  -  S(StP  Sj  .  .  .  4';*), 


where 


V  -  s„  (-  o/2  +  iit)-,  i  =  | ,  2,  3  . . .  2k. 

Let  us  study  the  statistical  properties  of  process  F^  *•  F(tm), 

m  ■  1,  2  ...  N.  Condition  (1)  indicates  that  the  telegraphic  signal, 

received  in  a  random  one-dimensional  background  cross  section,  is  a 
Markov  process.  Process  by  virtue  of  its  structure  resembles  the 

telegraphic  signal  and  differs  from  it  only  in  that  its  fronts  are 
smoothed  and  in  :hat  it  takes  both  positive  and  negative  values.  There¬ 
fore  we  shall  also  examine  process  as  a  random  Markov  sequence  with 

possible  f^  states  and  with  consequent  transition  probabilities.  From 
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state  fQ  the  only  possible  transition  is  to  states  fg,  f_^,  with  prob¬ 
abilities  a.' ,  oJ  2.  ail  respectively  (a'  *  1  -<X).  From  state  f^,  pro¬ 
cess  with  probability  ^  turns  into  and  with  probability  ■ 

1  -  p  remains  f^.  From  states  f  ^ ,  where  j  ^  1;  j  k;  f  4  2k -1,  with 
probability  1  a  transition  takes  place  to  state  fj+1*  Finally,  from 

state  £2^  i ,  the  process  with  probability  1  turns  into  f q.  Am  logout 
relationships  are  true  for  states  f_^,  where  j  >  0. 

Since  A  and  F  in  the  aggregate  form  a  Markov  process,  the  3-d:.- 
meni onal  vector  r 

«*,=0V,  A„  Fr) 

alsc  forms  a  Markov  process.  Therefore  the  detection  problem  obtained 
here  in  discrete  form  fits  totally  into  the  problem  scheme  used  in  the 
detection  of  a  weakly  determined  signal,  examined  in  reference  [5].  In 
Accordance  with  the  result  obtained  in  reference  [5],  the  optinal  re-  [6 
ceiver  consists  of  circuits  connected  in  series,  located  right  after 
the  init  to  eliminate  the  constant  component  of  Interference. 

1.  Of  a  nonlinear  transformation  effecting  "identification"  of  the 
inte  ference  and  transforming  input  realization  yr  into  the  sequence 


the  estimate  of  the  current  interference  value 

is  the  conditional  expection  value  ,  computed  on  the  basis  of  all 

observations  known  to  the  currei.t  time  moment  and  in  the  assumption 
that  the  signal  is  absent; 
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1  t 


l0(°i>  fj\y,) 


lJ(Ar~  a,:  F, 


are  a  posteriori  probabilities  determining  the  conditional  expectation 


EW,\y V)- 


2.  Of  a  linear  filter  consistent  with  the  form  of  the  signal  re¬ 
ceived  . 

3.  Of  a  threshold  unit. 

The  nonlinear  transformation  (6)  plays  the  major  role.  Since 
vector  (y  ,  A  ,  F  )  is  a  Markov  process,  the  a  posteriori  probabilities 
that  figure  in  (6^ 


»(«,;  fj  y,);  i  =  1.  -•••/>; 

J—  -(2*-  1)  ...  (2Ar  —  1) 


forn  a  conditional  Markov  process  [6].  This  very  important  property 
makes  it  possible  to  determine  the  conditional  probabilities  without 
the  need  of  remembering  all  the  current  values  of  vector  y  .  The 
method  of  computing  these  probabilities  has  been  quite  well  developed 
[6]  consequently,  omitting  intermediate  calculations,  we  shall  give 
the  final  relationships 


» («,;  /., \h  **  («'» (-«  I  /,-•)  ■: 

f  uk  .i)| /•*•-! I  >•/•-*)  1 ; 

®(«<;  Ail?,)  -- 

-  V  « (a- fjy,. ,),  (7) 


•(«*/;  Ail  jv)5® 


/-i 

h 


«(<v*  Aif-n1  >v-i; 


l>\  ,  I  Fk;  l-f-k-r  1. 

'  « («j;  A*  I  y,)  = - 7^-  -=—  !°>  l<h; 

A «.-!)  ly^-i)  +  ?'*  (<*,;  A*  I  >v  -O) ; 

»(«i:  A(*+uly,)— - 7, - w  < 

X(<v,  A*|yfi); 


where 


while  the  norming  factor  I  is  determined  ; rom  tl e  equation 

V  '2  *  («,;  AW—  '■  <S) 

As  seen  from  (8)  and  (9),  the  a  poste  iori  probabilities 

«fa;  //!?,) 

are  computed  recursively  from  the  current  .bservi'.tion  of  y  and  from 
the  values  of  these  same  probabilities  fou  td  for  the  preceding  time 
moment  r-1. 

Analysis  of  the  sensitivity  and  anti-  .nterft rence  capabilities  of 
the  detection  procedure  proposed  is  beyond  the  s<  ope  of  the  present 
article. 
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